QUANTUM CHAOS: A DECOHERENT DEFINITION. 
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ABSTRACT: We show that the rate of increase of von Neumann entropy computed 
from the reduced density matrix of an open quantum system is an excel·lent indicator 
of the dynamical behavior of its classical hamiltonian counterpart. In decohering quan- 
tum analogs of systems which exhibit classical hamiltonian chaos entropy production rate 
quickly tends to a constant which is given by the sum of the positive Lyapunov exponents, 
and falls off only as the system approaches equilibrium. By contrast, integrable systems 
tend to have entropy production rate which decreases as í" 1 well before equilibrium is 
attained. Thus, behavior of quantum systems in contact with the environment can be 
used as a test to determine the nature of their hamiltonian evolution. 



I. INTRODUCTION 



Ever since the inception of quantum theory, the issue of the correspondence between 
quantum and classical has been at the center of interest. Perhaps the best known (and most 
durable!) problem arises in the context of measurements 1 . The difficulty of representing 
the Universe as a whole - including us, the observers - by means of quantum theory has 
been known to the forefathers of quantum physics, and (in spite of the significant recent 
progress) continues to be hotly debated 2 '^. The problem of the correspondence between 
the quantum and the classical in quantum analogs of systems which classically exhibit 
dynamical chaos has come into focus within the past two dècades 4 '^. It is perhaps 
best illustrated by the fact that, in course of the hamiltonian evolution, quantum and 
classical versions of the same system begin to exhibit significant discrepàncies between the 
expectation vàlues of the same quantity on a relatively short timescale: 

t x oc À _1 log(xSp/h). (tchi) 

Above À is the Lyapunov exponent, x is the scale over which the potential becomes signif- 
icantly nonlinear, and Sp characterizes the initial spread of the wavepacket in the units of 
momentum. 

This timescale is only logarithmically dependent on the value of the Planck constant 
%. Thus, it is unconfortably short, even for macroscopic systems: If one were to take this 
prediction on the face value one would anticipate that in our Universe quantum chaotic 
systems should stop obeying classical laws after a few dynamical timescales (which is 
typically the order of the inverse Lyapunov exponent A _1 ). By contrast, quantized regular 
(integrable) systems tend to follow predictions of classical mechanics for much longer time 

- on a timescale of the order of where (3 is some (positive) characteristic power. This 
failure of the correspendence principle for chaotic systems has even led some to wonder 
whether quantum theory can be the fundamental theory of our Universe which - after all 

- seems to follow classical mechanics at the macroscopic level 8 . 

The purpose of this paper is to provide a brief summary of the relevant aspects of 
the process of decoherence - which was introduced 9 ' 10 to deal with the transition from 
quantum to classical in quantum measurements - and to show how it helps resolve the 
problem of quantum-classical correspondence in the context of chaos. As an important 
corollary of this discussion we will conjecture a simple method to diagnose chaos in the 
fully quantum system. We shall base it on the behavior of the von Neumann entropy 
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production rate of systems coupled to an environment. In cases where the classical system 
is chaotic, von Neumann entropy is conjectured to increase at a rate given by Lyapunov 
exponents in its decohering quantum analog. By contrast, in quantum analogs of regular 
systems entropy will grow at a rate which will asymptotically tend to zero well before the 
system reaches equilibrium. 

What is remarkable about this result 11 is that the entropy production rate - follow- 
ing the initial onset of decoherence which ocurrs on a decoherence timescale which is 
essentially independent of the system's self-hamiltonian, but dependent on the strength 
and nature of the coupling with the environment 1 ^ and on the form of the initial state^' 14 
- tends to be dictated by the dynamics of the system rather than by the type or the 
strength of the coupling. Independence from the strength and nature of the coupling holds 
for a wide range of parameters in spite of the fact that the ultimate cause of irreversibility 
is precisely the coupling with the environment. This behavior mirrors classical intuition 
about the nature of chaotic systems: Their evolution is - in contrast to regular systems - 
unpredictable. We show that this unpredictability carries over into the quantum domain, 
provided that systems which are operi are investigated. This result can be therefore re- 
garded as an additional indication that the correspondence between quantum and classical 
dynamics should be sought only with the assistance of environment - induced superse- 
lection, the consequence of the process of decoherence. Furthermore, entropy production 
rate in a decohering quantum system can be regarded as a diagnòstic tool: Rate of in- 
crease of entropy can distinguish chaotic and regular quantum evolutions, thus providing 
a completely quantum definition of quantum chaos. 



2. DECOHERENCE 



Decoherence and its relation with quantum measurement are not the main subjects 
of this paper. We shall summarize decoherence only very briefly with the eye on its 
significance to the subject of the transition from quantum to classical in the context of 
quantum chaos. A more complete review can be found elsewhere 1 ^' 1 ^ 

Decoherence is the process of loss of (phase) coherence by the system caused by the 
interaction with the external or internal degrees of freedom which cannot be followed by 
the observer and are summarily called 'the environment'. Different states in the Hilbert 
space of the system of interest show various degrees of susceptibility to decoherence. States 
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which are least susceptible (i.e., take longest to decohere) form the preferred basis (also 
known as the pointer basis, in the context of quantum measurement) 9 ' 15 ' 16 ' 17 . 

Preferred states are singled out by the interaction between the system and the environ- 
ment. In idealized discussions of quantum apparatus, complete immunity to decoherence 
can be guaranteed for the eigenstates of the pointer observable which commutes with the 
total hamiltonian (i.e., self hamiltonian plus the interaction hamiltonian) . Hence, pointer 
observable is conserved in spite of the interaction with the external degrees of freedom 9 . 

Monitoring by the environment is a useful way of thinking about the emergence of the 
preferred set of states and about the process of decoherence in general. It can be shown 
that the interaction with the environment can be regarded as a continuous measurement 
of the pointer observable. As a consequence, the environment acquires a record of pointer 
observable. Its states become correlated with the preferred pointer states: Quantum state 
of the complete (system plus environment) object can be written as; 

\ï& >= ^^ai\o~i > \t{ > , (totalwf) 

í 

where the states of the environment correlated to the eigenstates of the preferred observable 
become (to an excellent approximation) orthogonal (i.e., < q|ej >oc Sij) as a result of the 
interaction with the system. 

As the environment continuously acquires the record of the states {|c^ >}, other states 
(linear superpositions of the preferred states) are unstable. Quantum coherence of super- 
positions of {\ai >} is lost: When a system is prepared in a state; 

|$ >= ^^ai\ai >, (systemwf) 

i 

it will rapidly decay into a density matrix which is always diagonal in the same (preferred) 
basis: 

|$ >< $| — > p sa ^2 \ a i\ 2 \ a i >< a í\- (trànsit) 

i 

Hence, the system behaves as if an effective superselection rule precluding existence of 
superpositions between the eigenstates of the preferred basis was in place. 

Environment induced superselection rules effectively outlaw arbitrary superpositions. 
Thus, even though the superposition principle is vàlid in a closed quantum system, it is 
invalidated by decoherence for systems interacting with their environments. All of the 
macroscopic quantum systems we encounter in our everyday existence, as well as our own 
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memory and information processing hardware (e.g., neurons, etc) are macroscopic enough 
and sufficiently strongly coupled to the environment to be susceptible to decoherence, 
which will eliminate truly quantum superpositions on a very short timescale. This process 
is absolutely essential in the transition from quantum to classical in the context of quantum 
measurements (where the classical apparatus tends to be very macroscopic) although res- 
olutions based on decoherence may not be easily palatable to everyone (i.e., see comments 
on decoherence in the April 1993 issue of Physics Today). 

The timescale on which decoherence takes place can be estimated by solving a spe- 
cific example: a one dimensional particle moving in a potential V(x) coupled through its 
position with a thermal environment - e.g. with a collection of harmònic oscillators at a 
temperature T. Under the appropriate assumptions one can derive an equation for the 
reduced density matrix of the preferred particle. In the position representation it reads 18 : 



vori Neumann eq. relaxation decoherence 



P = ~[H ,p] ~ §^)P ~ -—^—^-yfp (mastereq) 



p—FORCE—'VV p=—jp classical phase space 

Above, we have indicated the role of the three terms which constitute the master 
equation in the so-called high temperature limit - that is, in the case when the thermal 
excitations of the environment dominate the effect of the environment, and the effect of 
the vacuum fluctuations can be neglected. 

Mathematically, classical limit is often associated with the size of Planck's constant. In 
the limit h — > the very last term of the above equation becomes dominant. To understand 
its effect, let us write an explicit solution of (mastereq) approximating the right hand side 
by the last term only. In that case 

p(x, y, t) = p(x, y, 0) exp(-(x - y) Dt/h ) (rhot) 

Above D = 2rri'jk B T. 

It is now apparent that the evolution under nothing but the decoherence term leaves 
the diagonal of the density matrix in the position representation essentialy uneffected. By 
contrast, for x ^ y the density matrix will decay exponentially in a decoherence timescale 12 ; 

r D = T- 1 ^^ = r R (^f , (decüme) 
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where À^j? = (fi 2 /2m/c^T) 1 / 2 is the thermal de Broglie wavelength and = 7 1 is the 
relaxation timescale. 

Two remarks are in order: (i) The decoherence timescale tjj is much shorter than 
the relaxation timescale tr for all macroscopic situations, as typical thermal de Broglie 
wavelengths of macroscopic bo dies are many orders of magnitude smaller than macroscopic 
separations Ax. (ii) The devastating effect of decoherence on superpositions of position 
can be traced back to the preferential monitoring of that observable (x) by the environ- 
ment, which was coupled to the position of the system of interest. This also tends to be 
the case in general: Interaction potentials depend on position and, therefore, allow the 
environment to monitor ^9,13,15,16 ^ g a resu i^ Q f the action of the decoherence term, the 
vast majority of states which could in principle describe the system of interest would be, in 
practice, eliminated by the resulting environment - induced superselection. Only localized 
states will be able to survive. They will form a preferred basis. For (even though they will 
be in general still somewhat unstable under the joint action of the self-hamiltonian and 
the environment) they will be much more stable than their coherent superpositions. This 
can be gauged by estimating the timescale characterizing the rate of entropy production. 
For the preferred states this timescale will be relatively long, determined by the dynamics 
and relaxation. For example, in an underdamped harmònic oscillator the preferred states 
turn out to be the familiar coherent states 14 : Oscillator dynamics rotates all of the states, 
which, in effect, translates spread in position into spread in momentum (and vice versa) 
every quarter period of the oscillation. As a result, coupling to position can be quite faith- 
fully represented in the "rotating wave approximation" which makes the master equation 
symmetric in x and p ^. Hence, coherent states will minimize entropy production. By 
contrast, for superpositions of coherent states entropy production will happen on a very 
much shorter decoherence timescale. 

More general dynamics tends to deform states in the Hilbert space of the system. 
Chaotic dynamics is especially effective in this, as it reflects exponential stretching and 
squeezing of phase space distributions. Thus, a regular patch in the phase space will be 
relatively quickly (on a Lyapunov timescale) deformed into something which will be no 
longer regular. By the same token, localized preferred states will tend to be stretched into 
non-local superpositions. Thus, the form of a typical state on the diagonal of the density 
matrix of a chaotic system will be a matter of compromise between the chaotic dynamics 
and decoherence. We shall discuss the nature of this compromise in the next section. 
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3. DECOHERENCE VS. EXPONENTIAL INSTABILITY 



Phase space provides a natural arena to study the consequences of the chaotic dynamics 
and its interplay with decoherence. Master equation (mastereq) can be translated into 
an equation for the Wigner distribution. The resulting equation consists of the Wigner 
transform of the commutator (which results in the so-called Moyal bracket) plus additional 
terms representing relaxation and decoherence: 

W = {H, W} MB + 2<yd pP W + Ddl p W (wignereq) 

where the first term on the right hand side denotes the Moyal bracket, which can be written 
in terms of the Poisson bracket as { , }mb = 2ism(h{ , }pp/2i)/h. 

We will be interested in the regime in which the coupling to the environment is suffi- 
ciently weak so that the damping (represented by the second term in (wignereq)) is neg- 
ligible. This is the so-called "reversible classical limit" 12 ' 15 ' 13 which in integrable systems 
yields reversible classical trajectòries but still eliminates non-local superpositions (this 
limit is achieved by letting 7 approach zero but keeping D constant so that decoherence 
continues to be effective). In this limit, and in the case where the potential is analytic, 
equation (wignereq) can be rewritten as 20 : 

W = {H, W} PB + Ç { fJl~^ 2 A 2n+1) H 2n+1) W + Dd 2 pp W. (wignereq2) 

Thus, Liouville flow in the phase space (and, therefore, classical dynamics) is obtained 
from the bàsic quantum picture as long as the corrections appearing in (wignereq2) are 
negligible. However, in a chaotic system evolution of the Wigner function generated by 
the Poisson bracket takes it quickly into the regime where Poisson bracket alone does 
not sufhce. This is because chaotic systems exhibit exponential sensitivity to initial con- 
ditions. Consequently, a phase space patch corresponding to a Wigner distribution will 
be exponentially stretched in the unstable directions corresponding to positive Lyapunov 
exponents. As the volume in phase space corresponding to W must be preserved, this 
will result in exponential shrinking in other directions. Consequently, derivatives of the 
Wigner function with respect to momentum (which enter into the correction term) will 
exponentially increase, so that after a time which is logarithmic in h these initially small 
terms will become comparable to the Poisson bracket and Liouville dynamics will cease 
to be an accurate approximation. This argument leads, in fact, to a demonstration of 
equation (tchi), as the reader is encouraged to verify. 
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One can regard this breakdown of the Liouville dynamics as a consequence of the loss 
of validity of a classical formula for the force in terms of the gradient of the potential 
V(x) (which is implemented in the Poisson bracket). As the Wigner function becomes 
more squeezed in momentum, by virtue of Heisenberg's uncertainty principle it spreads 
in position, and it begins to coherently sample increasingly large regions of the phase 
space. This process results in the domination of the evolution operator by the quantum 
forces when the extent of the wavefunction in space becomes comparable with the scale of 
nonlinearity, which for the various terms in equation (wignereq2) is given by: 

Xn = (d x V/d ( Í n+l) V) l l 2n . (chin) 

How can decoherence help reestablish the quantum - classical correspondence? Let us, for 
the moment, keep just the Poisson bracket and the diffusion term. Then, in the neigh- 
bourhood of any point, equation (wignereq2) can be easily expanded along the unstable 
(À+ > 0) and stable (\T < 0) directions in phase space (X^(\ r + = 0)- Diffusion 
will have little influence on the evolution of W along the unstable directions: After the 
possible initial (decoherence timescale) transient, W will be stretched simply as a result 
of the dynamics, so that the gradients along these directions will tend to decay anyway, 
without assistance from diffusion. By contrast, squeezing which occurs along the contract- 
ing directions will tend to be opposed by the diffusion. This will lead to a steady state 
with the solution asymptotically approaching a Gaussian with a half-width given by the 
critical dispersion: 

a c . = 2Di/\X i I (sigmacrit) 

where À~ is the (negative) Lyapunov exponent along the stable direction and is the 
diffusion coefficient along the same direction. Below, we will assume that the diffusion 
is isotropic (as would be the case in the rotating wave approximation) . Thus, after some 
time (and in the absence of folding - the other aspect of chaos which we will discuss below) 
the Wigner function will evolve into a multidimensional "hyper-pancake," still stretching 
along the unstable directions but with its width limited from below in the stable directions 
by equation (sigmacrit). 

At this stage, entropy will be approximated by the logarithm of the effective volume 
of the hyper-pancake. As its extent in the stable direction is fixed by the critical width 
(sigmacrit), its volume will tend to increase at a rate given by the positive exponents. 
Consequently, 

HttJ2 X t- ( hdot ) 

i 
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This constant rate will set in after a time larger than the decoherence timescale (for smaller 
times the entropy production can be even more ràpid) and after a time over which the 
initial Wigner distribution becomes squeezeed by the dynamics to the dimension of order 
of the critical dispersion a Ci . Equation (hdot) will be vàlid untill the pancake fills in the 
available phase space and the system reaches (approximately) uniform distribution ove the 
accessible part of the phase space, that is after a time defined by; 

t eq = (H eq /H )/H, (eqtime) 

where Hq is the initial entropy, and H eq is the entropy uniformized by the chaotic dynamics. 

Astute reader will note that H eq above need not be a true equilibrium entropy with the 
temperature given by T. Rather, it will correspond to dynamical quasi-equillibrium - the 
approximately uniform distribution over this part of the phase space which is accessible to 
the chaotic system as a result of its dynamics. 

The corresponding timescale will have a similar dependence on % as the timescale t% 
defined by Eq. (1). This is because entropy is approximately given by the logarithm of 
the volume of the phase space over which the probability distribution has spread in the 
units of Planck constant. Nevertheless, ty and t eq depend on rather different aspects of 
the initial and final state, and one can expect t x to be be typically a fraction of t eq . 

By contrast, in integrable systems stretching of the corresponding hyper-pancake in 
phase space will proceed only polynomially. Thus, even when it will get to the stage 
at which, in the contracting direction, diffusion will beco me important, stretching in the 
unstable direction will be only polynomial (rather than exponential) . Consequently, the 
volume of the hyper-pancake will increase only as some power of time. Hence, the entropy 
will grow only logarithmically as the entropy production rate will fall as H oc 1/t: It will 
take exponentially long to approach dynamical quasi-equilibrium. 

This difference in behavior between chaotic and integrable open quantum systems is 
striking and can be used as a defining feature of quantum chaos. 

4. QUANTUM-CLASSICAL 
CORRESPONDENCE IN CHAOTIC SYSTEMS. 

The failure of Ehrenfest theorem in chaotic systems is the consequence of the expo- 
nentially unstable Liouville flow which compresses Wigner function into an exponentially 
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narrowing pancake. As the momentum becomes progressively squeezed - which makes it 
less and less uncertain - the spatial extent of the coherent quantum wavefunction will ex- 
ponentially increase until it eventually simultaneously samples much of the potential well. 
By then the force is no longer given by a gradient of the potential: The wavefunction is too 
non-local for such a formula. It would not be even clear where (whithin the spatial sup- 
port of the wavefunction) one should compute such gradient. Instead, a more complicated 
formula, the Moyal bracket, is needed. 

Decoherence limits the extent over which the wavefunction can remain coherent. This 
is because a finite minimal dispersion in momentum (sigmacrit) corresponds to quantum 
coherence over distances no longer than: 

/ = h/o c = h(2D/X)- 1 / 2 . (Icrit) 

Thus, when the scale (chin) on which nonlinearities in the potential are significant is small 
compared to the extent of the wavefunction 

X <C / (chivsl) 

decoherence will have essentially no effect. Evolution will remain purely quantum and will 
be generated by the full Moyal bracket. 

By contrast, when the opposite is true, the evolution will never squeeze Wigner dis- 
tribution function enough for the full Moyal bracket to be relevant. Poisson bracket will 
suffice to approximate the flow of probability in phase space. The inequality character- 
izing this case can be written in a manner reminiscent of the Heisenberg indeterminacy 
principle: 

h <C x°~c- (ineqclass) 

That is, as long as decoherence keeps the state vector from becoming too narrow in mo- 
mentum, it will also prevent it from sampling the potential coherently over distances on 
which V(x) is noticeably different from linear. Hence, local gradients will suffice in the 
evaluation of the forces - Poisson bracket is all that is required. 

There is one more interesting regime where the chaotic motion is dynamically reversible 
(that is, H = 0) even if the system satisfies inequality (ineqclass). This happens when the 
initial patch in phase space is large (volume much larger than the Planck volume - initial 
entropy larger than a single bit) and regular. Then the initial stage of the evolution will 
proceed reversibly, in accord with the Poisson bracket generated flow. Decoherence will 
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have little effect. This is because its influence will set in only as the dimension of the 
Wigner distribution in the contracting dimension will approach the critical dispersion a c : 
In a simple example the entropy production will increase as: 

H = X- - I — — - (hdotapp) 

(l + - 1) exp(-2Aí)) 

So far, we have not taken into account (or, at least, not taken into account explicitelly) 
the other major characteristic of chaos: In addition to exponential instability, chaotic 
systems "fold" the phase space distribution. While this problem may require further study, 
we believe that the fundamentals of folding are already implicit in the above discussion: 
Folding will happen on the scale \ °f nonlinearities in the potential (which will typically 
- but not always - coincide with the size of the system, as it is defined by the range of 
its classical trajectory). Hence, preventing the system from maintaining coherence over 
distances of the order of x wm a ls° ascertain its classical behavior in course of folding. 
There will simply be no coherence left between the fragments of the wavepacket which 
will come into proximity as a result of folding, if they had to be separated by distances 
larger than / in the course of the preceding evolution. Thus, folding will proceed as if 
the system was classical, but with a proviso: After sufficiently many folds the distribution 
function (which in the stable direction cannot shrink to less than a c ) will simply fill in 
the available phase space. This will be achieved in the previously defined equilibrium 
timescale t eq . These conclusions are consistent with the numerical studies of quantum 
maps corresponding to open quantum systems such as the "standard map" carried out by 
Graham and his coworkers 21 . 



5. SUMMARY. 



We have used the entropy production rate in a decohering quantum system to char- 
acterize the nature of its evolution. Classical unpredictability - the essence of dynamical 
chaos - was shown to beget quantum unpredictability, quantified by the ràpid entropy pro- 
duction on the Lyapunov timescale. By contrast, dynamics of integrable systems leads to 
only gradual (polynomial) spread of the patch of the phase space corresponding to the state 
of the system. As a result, a much slower evolution towards dynamical quasi-equilibrium 
(and a relatively good predictability in spite of the coupling to the environmnet) charac- 
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terize quantum analogs of classically integrable systems. This distinction is conjectured to 
be a diagnòstic of the dynamical nature of the system. 
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Figure 1: The evolution of a "patch" in the phase space in a system with an exponential 
instability: 

(a) The case when the system is isolated and only exponential stretching in the unsta- 
ble direction as well as the corresponding shrinking in the complementary direction take 
place. The decrease of the dimension of the patch in momentum results (through the 
Heisenberg indeterminacy relation) in nonlocality, which leads to non-classical corrections 
to the expression for the force (resulting in the Moyal bracket). 

(b) When the system is open, decoherence prevents the dispersion of momentum from 
shrinking to less than te critical dispersion a c . Critical dispersion characterizes the steady 
state set by the competition between the dynamics (which attempts to narrow the patch, 
as it was shown in Fig. la) and decoherence, which is associated with the diffusion operator 
(which attempts to spread the patch). When decoherence is sufficiently effective, the spatial 
extent of the coherence of the wavepacket (given by / = h/a c , (lcrit)) will be sufficiently 
small so that the harmònic approximation to the potential will be accurate. And in such 
linear regime Moyal bracket and Poisson bracket coincide. Therefore, classical dynamics 
can be recovered even for chaotic systems. 



14 



